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ABSTRACT 

The size of the shock-layer governed by a conservation law is studied. The conservation 
law is a parabolic reaction-convection-diffusion equation with a small parameter multiplying 
the diffusion term and convex flux. Rigorous upper and lower bounding functions for the 
solution of the conservation law are established based on maximum-principle arguments. 
The bounding functions demonstrate that the size of the shock-layer is proportional to the 
parameter multiplying the diffusion term. 
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1. Introduction. Upper and lower bounding functions are presented. The bounds 
demonstrate that solutions to the singularly perturbed hyperbolic partial differential equa- 
tion 


(i) 


nr 1 & U 

: = It + 


d/M 

dx 


dx 2 


- R{u ) = 0 


with a pre-existing shock have shock-layers of width 0(c)- The analysis is performed in the 
style of Howes [4, 5, 6]. It begins with a multiple- scales asymptotic analysis. This provides 
the appropriate local scalings and indicates candidate forms for a bounding function. A 
bounding function is constructed from these candidates. Maximum principle arguments are 
then used to rigorously establish bounds for the solution. In this way, upper and lower 
bounds for solutions to equation (1) are established for any e. 

This presentation will concentrate on a comparison between solutions to (1) and solutions 
to the corresponding reduced equation 


( 2 ) 


PAU] := f + 


OW) 

dx 


R{U ) = 0 


obtained by setting e = 0. The result will be a bound on the difference between the solution 
to this reduced equation and the solution to equation (1). 

There are implications of this analysis for the computational aspects of the problem 
as well as for the physics modeled by conservation laws. The bounding functions result in 
an upper bound on the size of the shock-layer. They isolate the internal-layer region in 
which viscosity is important from the convection- dominated outer region. This validates the 
assumption that the lack of resolution of the physics in the shock-layer effects the solution 
in smooth regions very little under certain circumstances; thus, the hyperbolic equation (2) 
may be substituted for equation (1). 

This work sharpens the bounds of Howes [5, 7] which apply in a more general setting. 
The shock-birth region will not be studied here. 


2. Problem Specification. This paper is concerned with presenting bounds for the 
solution u to the quasilinear parabolic equation (1) on the domain 


( 3 ) 


D := {(s,t)|0 < x < b,0 < t < T}, 


subject to 

(4) u(x, 0) = g{x), 0 < x < 6; 

(5) u(0,t) — oc(t), 0 <t <T\ and 
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(6) u(b,t) = 0 <t<T. 

Let the portion of the boundary along which the data is specified be denoted by 


n := {(x, i)|0 < X < b, t = 0} |J{(x, f) |0 <t<T, x = 0,b}. 

The boundary and initial data are continuous and sufficiently smooth so that u is 
uniquely defined [1]. In particular, the boundary data a and (3 have derivatives bounded 
independent of e. Corner-layers in u are prohibited by assuming the compatibility conditions 

(7) oc(0) = ff(0), g(b)=/3( 0) 

(8) ^ + 3-m - R{g) = 0, for (*,«) = (0, 0); 
and 

(9) ^ - R(g ) = 0, for (x, t ) = (b, 0). 

For simplicity, it is assumed that all boundaries are inflow boundaries; hence, a(t)f'(a(t)) > 
o and mnm) < 0. Also assume that there is a single shock-layer in the initial data that 
is contained in an 0(eln(e)) neighborhood 7r 0 of (r 0 , 0), where T 0 is the location of the shock 
in U at t = 0. The domain of the initial viscous-layer is 7r 0 := {(x, t)\t = 0 and jx — r 0 | < 
eln(e)}. 

The solutions to the parabolic problem will be compared to a weak solution of the 
hyperbolic equation (2). Let U be the weak solution of (2) with boundary data (5-6) that is 
the solution to (1) in the limit as e > 0 tends to zero (denoted as e | 0). The initial condition 
will reflect a shock eminating from (f, x) = (OjTo). Thus, the initial condition for U is 

(10) U(x,0) = g 0 (x), 0 < x < b, 

where the difference g — go is zero except in 7r 0 . The relationship between g and go in 7To will 
be discussed in more detail in the proof of Theorem 4.4. 

Let the path of the shock in U be given by the curve (x,t) = (r(f),f). It is natural to 
describe the values of U at the shock as 

U R (t)= lim U(x,t) 

xir(t) 

and 


Ui ,{t) = lim U(x,t). 
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For ease of presentation we will assume that Ul > U R . The solution U to (2) will satisfy the 
entropy condition 

(11) U L (t) > S(t) > U R (t) 


where the speed S(t) o£ the shock is given by the Rankine-Hugoniot jump condition [8] 


( 12 ) 


_ mm - 

m ~ Vdt) - Mt) ' 


The entropy condition may be written as 

(13) MO = U U t ) - U R (t) > fio, 

for a constant fi 0 > 0 that is independent of e. It is assumed that the shock is part of the 
initial data and exists for the entire domain considered; thus, MO is defined for all i > 0. 
Notice that the initial condition for T may be obtained from g 0 and that T satisfies the 
ordinary differential equation dT(t)jdt = S(t). Thus, T is uniquely defined. 

3. Asymptotic Analysis. Some of the relevant physics is presented here. The discus- 
sion includes introducing the appropriate scales and equations when equation (1) is used to 
model a shock-layer. These results are obtained using a heuristic analysis that treats e as a 
small parameter. The results are exploited in the construction of the bounding function and 
are made more rigorous in Section 4. First some of the properties of the inviscid solution 

will be discussed. 

3.1. Derivation of Bounding Function. In this section we will derive a canonical 
form for the bounding function. The assumptions utilized in the derivation are used only to 
motivate the equation governing the canonical form. The assumptions are not necessarily 
valid. The validity of the bounding function will be established when it is used in Section 4. 

The canonical form for the solution of the viscous problem is derived here. Large gra- 
dients in the neighborhood of a shock-layer are resolved by using the spatial scale 1/e. This 
scale is combined with shock-following to obtain the internal-layer coordinates 


C = 



and r = t. 


The transformation defined by these coordinates is applied to equation (1) to obtain 


d f{u) _ Q2{l 

~er dc 



> 


where u(C, T ) = u{x,t). This suggests the regular expansion 


U — Uq + eui + «... 
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for u. We assume that this expansion is a priori valid in the shock layer and use identification 
in e to obtain the equation 




for Uq 


d( w d( d( 2 ' 

i. This equation is integrated with respect to ( once to obtain 

f(u 0 ) - S(t)uo - = const . 


To make this equation easier to solve, we approximate the term S — f juo by k( p , and we 
set the constant in the right-hand-side of the above equation to zero. The solution to this 
equation is the exponential 


(14) exp(— /c C P+1 )‘ 

Since the magnitude of the bounding function should be decaying away from the shock- 
layer, we assume k > 0, and substitute |£| for This is the candidate form for the bounding 
function. 

Guided by the form (14), let the general form for the bounding function uj be 

(15) u> = avt + 9e~ K * t 


Here x(r(i),x, e) is a linear measure of the distance between (x,t) and (T(t) y t). The power 
of the distance function is some positive p, and 6 is a non-negative function of t. 

4. Comparison Theorem. In this section we will show that the function -f U is an 
upper bound for u. The statement of this result is in the form of a comparison theorem. 
Several lemmas needed in the proof of Theorem 4.4 will be presented first. The first lemma is 
a maximum principle stated in the form most useful for the proof of the theorem. Lemma 4.2 
demonstrates how to choose the parameters in uj so that the maximum principle is satisfied 
within each of the regions fi 0 (*) := {x\0 <x< T(t)} and fl^i) := {x\ T(t) < x}. Finally, we 
see how the parameters in uj can be chosen to satisfy the maximum principle on (x, t) = (T, t). 

4.1. Maximum Principle. This lemma is a modification of the the Nagumo-Westphal 
Lemma [11] to include functions that are C°° except on a set of measure zero where they 
may be only C°. A condition on the spatial derivative replaces the condition involving a 
parabolic operator on this set. This lemma is a direct extension of the result by Nagumo 
and Westfal; thus, it is presented without proof. 

Lemma 4.1. Let z(x,t) be a continuous function that is differentiable except on a finite 
number of curves . Suppose 

(16) P[z) > 0 
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in the regions where z is differentiable, while z satisfies 


(17) 


(z)JX,t) = Urn 


dz(x,t) 

dx 


> lim 
X IX 


dz(x,t) 

dx 


M (*.0 


for curves ( x,t ) = ( X(t),t ) on which z is continuous but not differentiable. When these 
conditions arc satisfied along with 


(18) 


z > u for (x, t ) 6 II, 


then z >u throughout D . 

The implications of this lemma at discontinuities on the choice of the bounding function 
are demonstrated in Figure 1. The jump in the first derivative with respect to x is larger on 



Fig. 1. Discontinuous Upper Bounding Function 


the left of ( x,t ) = ( X,t ) than on the right. 

4.2. Analysis of u;: Continuous Regions. The results in this section require a con- 
vex flux function. Let Q(v) = • The assumption of convex flux means 

Q'{y ) > const 

for all v and for a positive constant independent of e where Q'(v ) = dC ^J v ^ • We also assume 
that Q'(v ) is bounded above. 

The shock speed S(t) and the values of Q(u) are related. The mean value theorem states 

S(t ) = Qifi) 


that 

(19) 
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for some intermediate function U(t ) in the open interval ]U R , U L [. We will assume the 
slightly more restrictive case of U in the closed interval [Ur + A, Ul - A], where A is a 
positive constant. This relationship between S and Q will be exploited in the proof of the 
following lemma. 


Lemma 4.2. There is an x \ close enough to T and an to small enough, such that for 
each region we may construct a specific form of uj to satisfy inequality (16). 

Proof. First some algebraic details are discussed. Let x — \x — T\/v. Differentiation of 
u) results in 


u t = av + 




e 


-K&P 


> 


and 






<*„ = ^~ 6p ( p . + e -Kf 

We use the mean value theorem to obtain f(U + 0 ;)* = Q(U + v)v x + (Q(U) + u>Q')U X) where 
Q' is evaluated at some function between U and U + u). The parabolic operator applied to 
the bounding function is 


P[U + u] = u t + Qu> x + QW X - eU xx - eu xx - R(U + u) 


or 

P[U + w] = 6 T\ + r 2 + —ZT 3 + — T4, 

V z V 

where 


and 


Ti = a — U xx + atQ'U x , 
t 2 = [6' + 6QU X }e~ K * r - R, 

r 3 = -6K 2 p 2 x 2( *-') ] e - K * 

r A = 6Kp[(S - Q )x p " 1 + (p - l)^- 2 ]e- Ksp . 
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We will consider two cases based on the two forms of cj to verify that (16) is satisfied: 
Case I. When i £ fi 0 we let u> — u> 0 , where 

ojq — act. 

We choose a large enough and T small enough so that T\ > 0. Inequality (16) holds in this 
case since T2, T3, 74 = 0. 

Case II. When x € fli we let u = uq, where 

uq = aet + pe~ Kt * . 

Here, 8 and k are positive constants independent of e, and p is a positive function of t. The 
term tj is positive from the choices in Case I; thus, it is sufficient to consider only T2, T3 and 
T4. We choose v » e. For example, we could use v = e|ln(e)|. Then, we assume that eo is 
small enough that With this choice of v, it is sufficient to show that r 4 is bounded below by 
a positive constant that is independent of e. 

First we will establish a lower bound for S{t) — Q(U + iu). These properties are based 
on the fluid dynamics properties of the problem. Namely, Q(u) is the speed at which the 
characteristics of u travel. The characteristics will be traveling faster than the shock for 
x to the left of T, and slower than the shock for x to the right of F . This means that 
S(t) — Q(U + u >) is positive for x > &o( K )- The relation between x and S — Q is depicted 
in Fig 2. (The precise shape of the curve cannot be known without more information about 
Q .) The location of &o(/c) moves closer to the origin as k increases. Also, multiplying S — Q 
by x reduces the variation near the origin. It is clear that we can choose k large enough such 
that 

1 + (5 — Q)x > 

□ 

We now have specific forms for the bounding function in each of the domains. Next we 
must show that we can satisfy inequality (17) and C° continuity with U + u> with the scaling 
and k of the previous lemma. 

4.3. Analysis of u: Discontinuities. First we will establish the C° continuity by 
choosing the parameters for the bounding function. Condition (18) will also be satisfied by 
the choices in this lemma. 

Lemma 4.3. We may choose the parameters in lji to satisfy (18) and such that u> is 

C°. 
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e 

Fig. 2. Example form for S-Q 

Proof. To obtain C° continuity let 0 = fi. 

The result will follow providing (18) is satisfied for (x,<) E II. This relation clearly holds 
except possibly for ( x,t ) € 7 r 0 . In this region it is necessary to impose some restriction on 
g — go. Assume that the viscous profile of the shock-layer is of an exponential type that is 
bounded by the function uj presented in this proof. This is not a severe restriction, since a 
layer with a profile like that of tanh[— (x/e)] (a solution to Burgers’ equation) is suitable. □ 

Theorem 4.4. Assume that g and go with their first and second derivatives are bounded 
independent of e except that go has a jump at (x, t ) = (To, 0) and g may have derivatives 
bounded depending on e for ( x,t ) E 7r 0 . Suppose u> is constructed fromu) 0 and u/j of Lemma 
4.2. Also assume that U x is continuous across (x,t) = (T, f). Then there is a positive eo such 
that 

(20) U + lu >u 

for ( x , t) E D when eo > e > 0. 

Proof. We will have only one case to consider based on imposing (17) on the curve 
(x,t) = (T(t),t). Consider a modified distance function x = |x — T + 8v\jv. Observe that 
for positive 6, the term ui x from equation (4.2) is negative and has magnitude 0(5). If we 
take 8 asymptotically close to zero, then we still satisfy this condition. Thus, by taking the 
limit, we may return to the original basis function since lim^ 0 ® = x. □ 

Remark 2. The constraint on U x across (x,t) = (T ,t) can be relaxed somewhat. This 
constraint was imposed so that 8 could be taken asymptotically close to zero in the proof. 
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With a more careful choice of 6 this constraint can be eased. 

Corollary 4.5. Under the conditions of the theorem, 

(21) \u-U\ = O (/xexp ~ ^ + 0(e), 

when eo > e > 0 and i/(e) is chosen as above. 

This result follows directly using symmetric arguments to obtain a lower bound. 

5. Implications. A direct result of this theorem is an upper bound on the size of the 
shock-layer. In this context, the shock-layer is defined as the region in which the solution to 
(1) differs from the solution to (2) by more than a specified amount. Namely, it is the region 
in which 

(22) \u~U\> a 

for some positive cr. As reflected in the following corollaries, there are different results 
depending on whether cr is a function of e or not. This result follows directly from the 
theorem. The following corollary extends this result to the case when a is independent of e 
and is a direct result of Corollary 4.5. 

Corollary 5.1. Suppose the conditions of Theorem 4. 5 obtain. Let (x(f),f) be the 
independent variables for which Inequality (22) is satisfied. If a is a constant independent of 
e then there is an €q small enough so that \u — U\ < <7 when 

|x(<) - r(()| < OM 

for e 0 > e> 0. 

When we define the shock-layer as the region such that 

|u(x,t) — U(x,t)\ > a, 
then there is a positive e 0 such that 

|x(t) - T(t)| < 0(eln(e)) 

for e 0 > e > 0. 

Physically motivated domain decomposition algorithms can be based on ideas presented 
in the analysis discussed herein. The computational domain can be partitioned into subdo- 
mains that have different physical behavior. Inside each of these regions different modeling 
equations (and hence different numerical methods) are used so computational effort is con- 
centrated on the relevent physics [10], [2]. Global error bounds have been developed for these 
domain decomposition methods [9]. In addition, the general idea of using different modeling 
equations can be extended to systems of equations [3]. 


9 


REFERENCES 


[1] J. R. Cannon, The One- Dimensional Heat Equation , vol. 23, Addison- Wesley Publishing Com- 

pany, Reading, Massachusetts, 1984. 

[2] Q. V. Dinh, R. Glowinski, J. Periaux, and G. Terrasson, On the coupling of viscous 

and inviscid models for incompressible fluid flows via domain decomposition, in First Inter- 
national Symposium on Domain Decomposition Methods for Partial Differential Equations, 
R. Glowinski, G. H. Golub, G. A. Meurant, and J. Periaux, eds., Philadelphia, 1988, SIAM, 
pp. 350-369. 

[3] Marc Garbey and Jeffrey S. Scroggs, Asymptotic-induced method for conservation laws , 

to appear in Proceedings for the Workshop on Asymptotic Analysis and Numerical Solution 
of Partial Differential Equations, H. Kaper and M. Garbey, eds., 1990. 

[4] F. A. Howes, Perturbed boundary value problems whose reduced solutions are nonsmooth , Indi- 

ana Univ. Math. J., 30 (1981), pp. 267-280. 

[5] , Multi- dimensional reaction- convection- diffusion equations , in Ordinary and Partial Differ- 

ential Equations, Proceedings of the Eighth Conference, Dundee, Scotland, B. D. Sleeman 
and R. J. Jarvis, eds., vol. 1151, New York, 1984, Springer- Verlag, pp. 217-223. 

[6] , Multi- dimensional initial- bound ary value problems with strong nonlinearities, Arch, for 

Rat. Mech. Anal., 91 (1986), pp. 153-168. 

[7] , Nonlinear initial-boundary value problems in Rn , Lectures in Applied Mathematics, 23 

(1986), pp. 259-273. 

[8] P. D. Lax, Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock 

Waves , SIAM, Philadelphia, 1973. 

[9] J. S. Scroggs, A physically motivated domain decomposition for singularly perturbed equations , 

SIAM Journal on Numerical Analysis, (to appear, 1990). 

[10] J. S. Scroggs and D. C. Sorensen, An asymptotic induced numerical method for the 

convection- diffusion-reaction equation , in Mathematics for Large Scale Computing, J. Diaz, 
ed., Marcel Decker, 1989, pp. 81-114. 

[11] W. Walter, Differential and Integral Inequalities , Springer- Verlag, New York, NY, 1970. 


10 




NASA Report Documentation Page 

Naionai AeftnavjKsHro ~ 

S;vkv ‘Vjr-i -cs'-iVO" 

1 Report No. 2. Government Accession No. 

NASA CR- 182080 

ICASE Report No. 90-52 

4. Title and Subtitle 

SHOCK-LAYER BOUNDS FOR A SINGULARLY PERTURBED 
EQUATION 


3. Recipient's Catalog No. 


5. Report Date 

August 1990 

6. Performing Organization Code 


7. Author(s) 


Jeffrey S. Scroggs 


8. Performing Organization Report No. 

90-52 


10. Work Unit No. 


9. Performing Organization Name and Address 

Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C , NASA Langley Research Center 
Hampton. VA 23665-5225 

12. Sponsoring Agency Name and Address 
National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23665-5225 


505-90-21-01 
11. Contract or Grant No. 

NAS1-18605 

13. Type of Report and Period Covered 

Contractor Report 

14. Sponsoring Agency Code 


15. Supplementary Notes 
Langley Technical Monitor: 
Richard W. Barnwell 


Submitted to SIAM Journal 
on Applied Mathematics 


Final Report — — 

16. Abstract 

The size of the shock-layer governed by a conservation law is studied. The 
conservation law is a parabolic reaction-convection-diffusion equation with a small 
parameter multiplying the diffusion term and convex flux. Rigorous upper and lower 
bounding functions for the solution of the conservation law are established based 
on maximum-principle arguments. The bounding functions demonstrate that the size 
of the shock-layer is proportional to the parameter multiplying the diffusion term. 


17. Key Words (Suggested by Author(s)} 


[ 18. Distribution Statement 


asymptotics, hyperbolic— parabolic partial 64 Numerical Analysis 
differential equations, maximum princi- 
ple 


Unclassified - Unlimited 


19. Security Classif . (of this report) 
Unclassified 


20 Security Classif. (of this page) 
Unclassified 


21 . No. of pages 
12 


| 22. Price 

A0 3 


NASA FORM 1626 OCT 86 


NASA* Langley, 1990 






